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TACTICAL  CONFLICT  DETECTION  AND  RESOLUTION  IN  A  3-D  AIRSPACE 


GILLES  DOWEKt,  CESAR  MUNOZ*,  AND  ALFONS  GESER§ 

Abstract.  This  paper  presents  an  algorithm  for  detection  and  resolution  of  air  traffic  conflicts  in  a  3- 
dimensional  (3-D)  airspace  for  two  aircraft,  namely  ownship  and  intruder.  A  conflict  is  a  projected  incursion 
of  the  intruder  aircraft  within  the  protected  zone  of  the  ownship.  A  solution  is  a  single  maneuver,  to  be 
performed  by  the  ownship,  that  effectively  keeps  the  required  minimum  separation  without  cooperation  of 
the  intruder  aircraft.  The  input  to  the  algorithm  is  the  state  information,  i.e. ,  horizontal  position,  altitude, 
ground  track,  and  vertical  and  ground  speed,  of  both  aircraft.  The  algorithm  outputs  a  set  of  solutions. 
Each  solution  modifies  only  one  state  parameter  of  the  ownship:  ground  track,  ground  speed,  or  vertical 
speed.  The  proposed  algorithm  is  suitable  for  formal  verification. 

Key  words,  conflict  detection,  conflict  resolution,  CD&R  algorithm,  3-dimensional,  collision  avoidance 

Subject  classification.  Computer  Science 

1.  Introduction.  One  of  the  main  elements  of  the  Free-Flight  concept  [11]  is  the  redistribution  of 
responsibilities  for  air  traffic  separation.  Under  Free-Flight  rules,  each  aircraft  with  an  appropriate  leve  of 
equipment  is  responsible  to  assure  separation  with  other  aircraft  in  the  vicinity.  To  support  this  mode  of 
operation,  several  automated  decision  support  systems  are  being  proposed.  In  this  context,  Conflict  Detection 
and  Resolution  algorithms  (CD&R)  are  designed  to  warn  pilots  about  an  imminent  loss  of  separation,  and 
to  assist  them  in  a  corrective  maneuver. 

In  this  paper,  we  present  a  tactical  CD&R  algorithm  in  a  3-D  space  for  two  aircraft.  In  CD&R-related 
literature,  tactical  algorithms  use  only  state  information  to  project  aircraft  trajectories.  Due  to  this  intention¬ 
ally  limited  source  of  information,  they  are  intended  to  be  used  with  short  lookahead  times  (a  few  minutes, 
typically  5-10)  during  which  aircraft  are  supposed  to  follow  straight  flight  paths.  Strategic  approaches,  in 
contrast,  use  intent  information  such  as  flight  plans,  and  uncertainties  such  as  weather  conditions.  They 
may  have  lookahead  windows  of  several  minutes  and  even  hours.  For  a  survey  on  CD&R  methods  see  [7]. 

The  input  to  our  algorithm  is  the  state  information,  i.e.,  3-D  position,  ground  track,  vertical  speed,  and 
ground  speed,  of  two  aircraft.  We  distinguish  one  aircraft  as  the  ownship  and  the  other  as  the  intruder.  Loss 
of  separation  between  the  aircraft  is  predicted  via  linear  projections  on  time  of  the  state  parameters.  In  case 
of  a  predicted  conflict,  the  algorithm  proposes  several  solutions  for  the  ownship.  Every  solution  is  a  single 
maneuver  that  effectively  avoids  the  conflict.  The  maneuver  modifies  only  one  parameter  of  the  ownship. 
This  constraint  produces  finitely  many  solutions,  simplifies  the  calculations  performed  by  the  algorithm,  and 
is  simple  to  conceive  and  to  perform  by  the  crew. 

In  Section  2,  we  briefly  survey  some  methods  for  conflict  detection  and  resolution  that  inspired  our 
approach.  In  Section  3,  we  present  our  approach  and  its  theoretical  support.  In  Section  4,  we  give  an 
algorithm  to  find  convenient  3-D  conflict  solutions.  A  prototype  implementation  is  described  in  Section  5. 
The  last  section  summarizes  our  work  and  suggests  lines  of  research  for  future  work. 
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2.  Conflict  Detection  and  Resolution  Supporting  Free  Flight.  Distributed  Air/Ground  Traffic 
Management  (DAG-TM)  [1]  is  a  set  of  conceptual  elements,  developed  within  the  Advanced  Air  Transporta¬ 
tion  Technologies  project  at  NASA,  that  defines  modes  of  operation  supporting  the  Free-Flight  concept. 

Safety  assessment  of  new  air  traffic  management  systems  is  a  main  issue  in  DAG-TM.  Prototype  tools 
such  as  the  Autonomous  Operations  Planner  (AOP)  are  being  developed  at  NASA  Langley  to  study  the 
feasibility  of  self  separation.  Systems  with  similar  goals  have  been  proposed  in  other  research  laboratories, 
e.g.,  the  Future  ATM  Concepts  Evaluation  Tool  (FACET)  [3]  at  NASA  Ames1  and  the  Airborne  Separation 
Assurance  System  (ASAS)  [6]  at  the  National  Aerospace  Laboratory  (NLR)  in  the  Netherlands.  All  these 
tools  implement  CD&R  algorithms. 

Standard  safety  assessment  techniques  such  as  testing  and  simulation,  although  useful,  have  serious 
limitations  in  new  systems  which  are  significantly  more  autonomous  than  the  older  ones.  Given  the  critical 
nature  of  the  problem,  we  believe  that  safety  statements  should  be  made  and  verified  formally,  and  that 
proofs  should  be  checked  by  machine. 

With  the  above  premise  in  mind,  we  have  examined  at  the  AOP  conflict  detection  and  resolution  proce¬ 
dures.  The  AOP  will  eventually  provide  both  tactical  and  strategic  CD&R  resolution.  In  what  follows,  we 
consider  the  conflict  detection  part  of  the  CD&R  AOP  algorithm,  which  is  an  adaptation  of  a  deterministic 
procedure  implemented  within  the  ground  NASA’s  Center/TRANCON  Automation  System  (CTAS)  [12]. 2 

2.1.  Conflicts  and  Protected  Zone.  Two  aircraft  are  said  to  be  in  conflict  if  their  vertical  separation 
is  (strictly)  less  than  H,  H  >  0,  and  their  horizontal  separation  is  (strictly)  less  than  D,  D  >  0.  A  body 
in  the  3-D  space,  called  protected  area ,  is  assigned  to  each  aircraft  such  that  a  conflict  is  equivalent  to  an 
intrusion  of  another  airplane  into  its  protected  area.  The  protected  area  forms  a  cylinder  (hockey-puck, 
nickel,  pizza)  of  altitude  H  and  radius  D  around  the  position  of  the  aircraft.  The  values  H  =  1000  ft  (feet) 
and  D  —  5  nm  (nautical  miles)  are  commonly  used. 

Note  that  the  boundaries  are  not  considered  part  of  the  protected  area.  We  will  see  later  that  this  choice 
enables  optimal  ownship  maneuvers  that  touch  the  boundaries  of  the  intruder  protected  area. 

2.2.  AOP  Conflict  Detection.  The  input  to  the  strategic  AOP  conflict  detection  algorithm  [8]  is  a 
set  of  trajectories,  one  of  which  is  the  ownship  trajectory.  A  trajectory  is  a  list  of  points,  called  nodal  points , 
which  are  assumed  to  be  joined  by  linear  segments.  Each  nodal  point  contains  the  intended  aircraft  state 
(position  and  velocity  vector)  at  a  given  time.  Since  conflict  detection  is  assumed  to  be  asynchronous,  the 
first  step  of  the  algorithm  is  to  synchronize  all  the  trajectories.  This  is  done  by  taking  time  steps  of  duration 
A  and  then  measuring  the  distance  between  the  trajectories  at  every  step  during  a  lookahead  period  of 
time.  If  after  n  time  steps  there  is  a  violation  of  the  ownship  protected  zone,  then  a  conflict  is  detected. 
We  illustrate  the  situation  in  Figure  2.1,  where  a  loss  of  separation  occurs  at  time  nA.  The  algorithm 
implements  several  heuristics  to  avoid  unnecessary  calculations. 

Since  the  algorithm  does  not  compute  the  actual  time  when  the  first  loss  of  separation  occurs,  the  choice 
of  A  is  crucial  in  this  approach.  Indeed,  if  A  is  too  large,  near  misses  can  occur  without  being  detected. 
CTAS  uses  A  =  10  seconds. 

The  synchronization  step  introduces,  in  fact,  a  time  and  space  discretization.  In  recent  work,  we  have 
discovered  that  discretization  makes  formal  verification  more  difficult  [4].  Moreover,  discretization  of  tra¬ 
jectories  can  lead  to  an  accumulation  of  modeling  inaccuracies  that  lead  to  imprecise  conclusions.  We 

1  FACET  is  a  CD&R  analysis  tool  rather  than  a  flight  deck  decision  support  tool. 

2 The  AOP  detection  algorithm  also  detects  conflicts  with  hazard  areas.  That  kind  of  detection  is  outside  the  scope  of  this 
paper. 
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Fig.  2.1.  Strategic  AOP  conflict  detection 

have  successfully  verified  a  conflict  alerting  algorithm  using  a  continuous  trajectory  model  [9].  Continuous 
trajectory  models  are  well  suited  to  formal  verification. 

The  approaches  exhibited  in  the  remainder  of  this  section  use  a  2-D  geometry. 

2.3.  Conflict  Detection.  Let  aQ,  vQ  and  vl  be  the  position  and  ground  speed  vector  of  the  ownship 
and  intruder  aircraft,  respectively,  at  time  0.  We  assume  that  the  ground  speeds  are  constants.  Separation 
is  lost  at  time  t  if  and  only  if  the  projected  distance  between  both  aircraft  at  time  t  is  strictly  less  than  D, 
i.e., 


\(a0  ~  ai)  +  t(v0  -  Vi)\2  <  D2  .  (2.1) 

This  constraint  has  solutions  for  t  if  and  only  if  Equation  2.2 

\(a0  -  ai)  +  t(v0  -  vl)\2  =  D2  (2.2) 

has  two  solutions  t\  and  1 2  such  that  t\  7^  £2-  The  two  solutions  correspond  to  the  times  when  the  loss  of 
separation  starts  and  ends.  Note  that  by  definition,  if  t\  —  1 2  no  lost  of  separation  occurs. 

The  above  procedure  is  used  in  both  ASAS  (modified  potential  field)  and  FACET  (Bilimoria’s  geometric 
optimization)  algorithms  for  conflict  detection  in  2-D  [6,  2]. 

2.4.  Modified  Potential  Field  Resolution.  If  a  conflict  is  predicted,  the  modified  potential  field 
approach  (originally  due  to  [5]  and  implemented  in  ASAS)  computes  the  time  of  closest  separation  r  between 
the  ownship  and  the  aircraft  (r  =  {t\  +  ^2) /2) .  Then,  a  new  speed  vector  for  the  ownship  is  calculated  such 
that  the  distance  at  time  r  between  the  aircraft  is  exactly  D.  The  solution  is  illustrated  in  Figure  2.2  where 
c?0(t)  and  afir)  are  the  ownship  and  intruder  projected  positions,  respectively,  at  time  r,  v’0  is  the  ownship’s 
new  speed  vector,  and  a'0(r )  is  the  new  projected  position  of  the  ownship  at  time  r. 

The  modified  potential  field  approach  does  not  solve  conflicts.  If  the  ownship  maneuvers  towards  the  new 
speed  vector,  and  no  further  action  is  taken,  then  there  will  be  a  conflict  after  time  r.  Without  cooperation 
of  the  intruder,  the  ownship  will  have  to  pursue  repeated  maneuvers  to  only  approximate  a  solution.  To 
completely  solve  the  conflict,  even  an  infinite  number  of  maneuvers  are  necessary. 

2.5.  Geometric  Optimization  Resolution.  In  this  approach,  proposed  by  Bilimoria  [2]  and  imple¬ 
mented  in  FACET,  the  intruder  is  considered  fixed  in  space,  and  the  ownship  position  a  and  velocity  vector 
v  are  taken  relative  to  the  intruder  state,  i.e.,  a  =  aQ  —  cti  and  v  =  vQ  —  A  new  relative  speed  vector  for 
the  ownship  solves  the  conflict  if  it  does  not  intersect  the  interior  of  the  intruder  protected  area.  Among 
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Fig.  2.2.  Modified  Potential  Field  resolution 


the  infinitely  many  new  speed  vectors  that  solve  the  conflict,  Bilimoria  chooses  those  which  minimize  their 
angle  to  the  original  speed  vector.  Such  speed  vectors  are  called  optimal 

As  illustrated  in  Figure  2.3,  any  optimal  solution  is  tangential  to  the  intruder’s  protected  zone.  Any 
other  solution  requires  a  greater  change  of  the  ownship  relative  ground  track.  Each  touch  point  A  and  B 
determines  a  new  ownship  relative  ground  track  which  is  optimal  under  certain  constraints.  For  instance  the 
target  point  A  is  optimal  under  the  constraint  that  only  a  left  turn  may  be  made.  The  length  of  the  speed 
vector  may  be  arbitrarily  chosen.  A  minimal  solution  is  proposed  in  [2]. 


Fig.  2.3.  Geometric  Optimization  resolution 


3.  An  Approach  to  3-D  CD&R.  The  Geometric  Optimization  algorithm  in  [2]  uses  a  2-D  geometry, 
i.e.,  it  detects  and  solves  conflicts  in  the  horizontal  plane.  In  the  Modified  Potential  Field  algorithm  imple¬ 
mented  in  AS  AS  [6],  3-D  conflicts  are  decomposed  into  horizontal  and  vertical  conflicts  which  are  detected 
and  solved  independently.  Then,  the  solutions  are  composed  to  obtain  a  3-D  maneuver.  This  approach  is 
appealing  for  its  simplicity,  but  it  is  rather  difficult  to  prove  correct. 

We  pursue  a  true  3-D  geometric  analysis  to  conflict  detection  and  resolution.  In  this  section  we  extend 
Bilimoria’s  horizontal  CD&R  approach  to  three  dimensions. 

Given  the  positions  and  speed  vectors  of  two  aircraft  (ownship  and  intruder),  we  compute  the  relative 
position  a  —  v0  —  Vi  and  relative  speed  v  =  v*0  —  vl  of  the  ownship  with  respect  to  intruder  referential.  We 
take  a  coordinate  system  where  the  origin  is  at  the  intruder  position  and  ax  <  0,ay  =  0  (see  Figure  3.1). 
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Fig.  3.1.  3-D  coordinate  system  and  protected  zone 

3.1.  3-D  Conflict  Detection.  The  relative  ownship  trajectory  L  is  the  half  line 

L  =  {a  +  tv  |  t  >  0}.  (3.1) 

We  say  that  there  is  a  conflict  between  the  ownship  and  the  intruder  aircraft  if  and  only  if  a  E  P,  where  P 
is  the  protected  zone 


P  =  {(x,  y,z)  \x2  +y2  <D2  and  -H  <z<  H}.  (3.2) 

The  two  aircraft  are  predicted  to  be  in  conflict  if  the  relative  trajectory  L  intersects  the  protected  zone  P. 
To  compute  the  time  interval  where  the  projected  intrusion  occurs,  we  compute  the  intersection  between  L 
and  the  boundary  of  P.  The  boundary  consists  of  two  parts:  (1)  the  lateral  surface  around  the  cylinder 

Pi  =  {(x,  y,  z)  |  x2  +  y2  =  D2  and  -H  <z<  H}  (3.3) 

and  (2)  the  top  and  bottom  bases 

P2  =  {(x,y,z)  |  x2  +  y2  <  D2  and  \z\  =  H}.  (3.4) 

Observe  that  we  count  the  top  and  bottom  circles  to  the  lateral  surface  and  not  to  the  bases. 

If  the  half  line  L  does  not  intersect  Pi  nor  P2,  then  the  two  aircraft  are  predicted  to  be  in  conflict. 
Otherwise,  we  compute  times  t\  and  t2  such  that  aircraft  will  be  in  conflict  during  the  interval  t\  <  t  <  t2. 

Assume  that  the  ownship  and  the  intruder  are  not  in  conflict.  The  following  algorithm  returns  true  if 
the  two  aircraft  are  predicted  to  be  in  conflict,  and  false  otherwise.  It  has  two  parameters:  a  =  (ax,0,az) 
and  v  =  (vx,vy,vz). 

1.  Case  vz  =  0.  This  case  corresponds  to  2-D  resolution.  Return  true  if  \az\  <  H,  vx  >  0,  and  the 
discriminant  D2{v2  +  v2)  —  a2v2  is  positive.  Otherwise,  return  false. 

2.  Case  vz  ^  0.  Let  t\  —  H~az  and  t2  =  ~H~az  •  These  are  the  times  when  the  trajectory  reaches 
the  altitudes  ±H.  Let  di,d2  denote  the  squares  of  the  horizontal  distances  to  the  intruder  at  times 
ti,t2,  respectively.  We  have  three  cases: 

(a)  Case  d\  <  D 2  and  d2  <  D2 .  There  are  two  intersections  with  the  bases  P2,  provided  that  the 
times  ti,t2  are  in  the  future.  Return  true  if  min(ti,t2)  >  0,  and  return  false  otherwise. 
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(b)  Case  d\  >  D 2  and  d 2  >  D2 .  In  this  case  there  are  no  intersections  with  the  bases,  but  there  may 
be  intersections  with  the  lateral  surface  Pi.  Return  true  if  vx  >  0  and  D2{v2  +  v 2)  —  a2v2  >  0. 
Otherwise,  return  false. 

(c)  Case  d\  <  D 2  and  d^  >  D2 .  We  have  one  intersection  with  a  base  at  time  t\.  Return  true  if 
t\  >  0.  Otherwise,  return  false. 

(d)  Case  d\  >  D 2  and  d^  <  D2 .  We  have  one  intersection  with  a  base  at  time  1 2.  Return  true  if 
£2  >  0.  Otherwise,  return  false. 

3.2.  3-D  Conflict  Resolution:  The  Target  Set.  We  assume  that  the  ownship  and  the  intruder  are 
not  yet  in  conflict,  but  that  they  are  predicted  to  be  in  a  conflict.  We  want  to  modify  the  relative  speed 
vector  v  to  a  vector  v'  such  that  the  half  line  L'  =  {a  +  v’t  \  t  >  0},  does  not  intersect  the  protected  zone. 

Among  the  various  solutions,  we  focus  on  those  that  modify  the  vector  vf  in  an  optimal  way,  i.e. ,  such 
that  L'  touches  the  boundary  of  P.  We  will  also  consider  the  special  solution  vf  =  0  (i.e.,  the  two  aircraft 
have  exactly  the  same  speed)  as  optimal.  Positive  multiples  of  optimal  solutions  are  again  optimal  solutions. 
Hence,  we  shall  first  characterize  the  directions  of  optimal  solutions. 

In  what  follows  we  assume  that  a  is  not  on  the  boundary  of  the  protected  area  P  (the  case  where  a  is 
on  this  boundary  must  be  handled  separately),  direction  of  an  optimal  solution  is  determined  by  one  target 
point  of  the  half  line  L'  with  the  boundary  of  P.  We  must  indeed  be  aware  of  half  lines  touching  one  of  the 
boundary  circles  of  the  lateral  surface  in  two  points  if  az  =  ± H.  Let  a  target  set  of  a  be  a  set  of  target  points 
such  that  for  every  optimal  solution  one  of  its  target  points  is  in  the  set.  The  kernel  of  our  3-D  conflict 
resolution  algorithm  is  the  computation  of  a  suitable  target  set. 

A  target  point  may  be  either  on  the  lateral  surface  Pi  or  on  the  bases  P2.  If  we  have  a  target  point 
on  one  of  the  bases,  then  V  must  touch  an  open  line  segment  of  the  base,  and  moreover  it  must  touch  the 
lateral  surface  at  two  points.  In  this  case  we  decide  to  take  one  of  the  latter  points  as  target  points.  So  all 
target  points  are  on  the  lateral  surface. 

A  target  point  satisfies  two  things:  (1)  it  is  on  the  boundary  of  the  protected  zone,  and  (2)  the  half  line 
V  passing  on  this  point  must  not  intersect  the  protected  zone  (only  its  boundary).  If  the  target  point  has 
coordinates  (x,  y ,  z)  the  first  condition  rephrases  x2  +  y2  =  D 2  and  —H  <  z  <  H.  The  second  is  that  for  all 
time  t  >  0,  (ax  +  tvx ,  tvy ,  az  +  tvz )  0  P. 

Our  main  goal  on  this  computation  is  to  remove  the  phrase  “for  all  time  t  >  0”  to  get  an  algebraic 
characterization  of  the  target  set.  Let  to  be  the  time  when  the  half  line  L'  crosses  the  target  point  (x,y,z), 
i.e.,  x  =  ax  +  v'xto,  y  =  v'yto ,  2  =  az  +  v'zto-  We  define  T  —  tjto  as  the  normalized  time  on  this  half  line.  The 
half  line  V  can  be  written  as  the  set  of  points  (ax  +  T (x  —  ax) ,Ty ,  az  +  T (z  —  az))  for  T  >  0.  The  condition 
that  this  half  line  does  not  intersect  the  protected  zone  is  rephrased:  For  all  T  >  0  either 

(ax+T(x-ax)f  +T2y2  >D2  (3.5) 

or 

\az+T(z-az)\  >H.  (3.6) 

Using  x2  +  y2  =  D2 .  Formula  3.5  yields 

(T  -  1)[(T  -  l)((s  -  ax)2  +  y2)  +  2 (D2  -  axx )]  >  0  (3.7) 

and  Formula  3.6  yields 

|*  +  (T-l)(*-a*)|  >H.  (3.8) 
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First  we  consider  the  target  points  such  that  —H  <  z  <  H  and  then  those  such  that  \z\  =  H. 

•  Case  —H  <  z  <  H.  For  some  time  T  in  an  interval  around  1,  Formula  3.8  does  not  hold,  i.e., 

\z  +  (T  -  l)(z  -  az)\  <H. 

Therefore,  Formula  3.7  must  hold  for  an  interval  around  1.  Thus,  we  must  have  D 2 3 4 5  —  axx  =  0  and 
hence 

x  =  D2 /ax. 

As  x2  +  y2  =  D2,  we  get 

y  =  ey 'D2  -  WF*?  =  sDVrt  -  D*/ax, 

where  e  is  ±1. 

As  the  target  set  must  be  on  the  boundary  of  the  protected  area,  we  must  have  — ax  >  D. 

It  is  easy  to  check  that  all  points  such  that  x  =  D2 /ax,  y  =  —eD^a\  —  D2 /ax  and  —H  <  z  <  H 
are  target  points. 

•  Case  \z\  =  H.  If  z  and  z  —  az  have  the  same  sign,  Formula  3.8  is  equivalent  to  T  >  1,  and  we  must 
have  for  all  T  <  1 

(T  -  1)[(T  -  l)((a:  -  ax )2  +  y 2)  +  2 (D2  -  axx)  >  0 

and  hence  D 2  —  axx  <  0,  i.e.,  x  <  D2 /ax.  Symmetrically,  if  2  and  z  —  az  have  different  signs,  we 
must  have  x  >  D2 /ax. 

The  relative  signs  of  2  and  z  —  az  depends  of  2,  which  can  be  ±i7,  and  the  position  of  az  with 
respect  to  ±H. 

We  analyze  all  the  cases. 

1.  Case  —H  <  az  <  H  and  —  ax  >  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,  y ,  z)  of  Pi  such 
that  x  —  D2  jax  or  (z  =  H  and  x  <  D2 /ax)  or  (z  =  —H  and  x  <  D2 /ax).  See  Figure  3.2. 


Fig.  3.2.  Target  set.  Case  —H  <  az  <  H  and  —ax  >  D 

2.  Case  az  <  —H  and  —  ax  >  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,  y ,  z)  of  P\  such  that 
x  =  D2 /ax  or  (z  =  H  and  x  <  D2 /ax)  or  (z  =  —H  and  x  >  D2 /ax).  See  Figure  3.3. 

3.  Case  H  <  az  and  —  ax  >  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,y,z)  of  P\  such  that 
x  =  D2 /ax  or  (z  =  H  and  x  >  D2 /ax)  or  (z  =  —H  and  x  <  D2 /ax).  See  Figure  3.4. 

4.  Case  az  <  —H  and  —  ax  <  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,  y ,  z)  of  P±  such  that 
2  =  —H.  See  Figure  3.5. 

5.  Case  H  <  az  and  —  ax  <  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,y,z)  of  P\  such  that 
z  —  PL.  See  Figure  3.6. 
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Fig.  3.3.  Target  set.  Case  az  <  —H  and  —ax  >  D 


Fig.  3.4.  Target  set.  Case  H  <  az  and  —ax  >  D 


6.  Case  az  =  —H  and  —  ax  >  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,  y ,  z)  of  P\  such  that 
x  =  D2 /ax  or  (z  =  H  and  x  <  D2 jax)  or  2:  =  —H.  See  Figure  3.7.  Notice  that  we  get  the  same  set 
of  directions  as  the  target  set  in  Figure  3.3.  These  directions  are  already  included  in  case  (2). 

7.  Case  az  —  H  and  —  ax  >  D.  In  this  case,  the  target  set  is  the  set  of  points  (x,y,z)  of  P\  such  that 
x  —  D2 /ax  or  2  =  H  or  (z  =  —H  and  x  <  D2 /ax).  See  Figure  3.8.  We  get  the  same  set  of  directions 
as  the  target  set  in  Figure  3.4.  These  directions  are  already  included  in  case  (3). 

4.  Constrained  Solutions.  In  the  previous  section,  we  got  an  infinite  set  of  target  points  (which 
defines  optimal  directions  for  vf).  A  new  ownship  speed  vector  v'0  can  be  calculated  from  any  point  in  the 
target  set  and  any  length  of  the  relative  speed  vector  v' .  Among  all  these  vectors,  representing  maneuvers, 
some  are  more  convenient  than  others.  In  this  section,  we  select  solutions  where  only  one  parameter  of  the 
(absolute)  ownship  speed  vector  is  modified,  i.e.,  ground  speed,  ground  track,  or  vertical  speed. 

As  we  have  seen,  the  target  set  is  a  subset  of  the  set  of  points  in  the  set  P1  =  P[  U  P2  (see  Figure  4.1), 
where 

P[  —  {(x,  y ,  z)  |  x2  Py2  —  D 2  and  2  =  eH  and  e  —  ±1}  (points  in  circles), 

P2  —  {(x,y,z)  |  x  =  D2 /ax  and  y  =  —sD yj —  D2 /ax  and  e  —  ±1}  (points  in  lines). 

To  get  constrained  solutions,  we  first  compute  the  solutions  on  this  set  satisfying  a  given  constraint.  For 
that,  we  assume  that  at  some  time  £,  the  point  (x,y,z)  =  (ax  +  t(vfox  —  Vix),t{v'oy  =  Viy),az  +  t(v'oz  —  Viz)) 
is  in  P1  and  we  proceed  by  case  analysis  on  P[  and  P2.  A  special  case,  when  v'  =  0,  is  also  considered.  We 
ignore  points  where  t  >  0.  Once  the  solutions  are  found,  we  check  whether  they  belong  to  the  target  set  or 
not. 


4.1.  Ground  Speed  Change  Only.  We  have  v'0  =  (k  vOXlk  voy,voz )  for  some  k  >  0.  We  must 
determine  the  possible  k  positive  such  that  at  some  time  £,  the  point  (x,y,z)  =  (ax  +  t(v'ox  —  Vix),t(v'oy  — 
Viy),az  +  t(v’oz  —  Viz))  is  on  the  target  set. 


8 


Fig.  3.5.  Target  set.  Case  az  <  —H  and  —ax  <  D 


Fig.  3.6.  Target  set.  Case  H  <  az  and  —ax  <  D 


1.  Points  on  the  circles  P[.  We  have 


2  i  2  7 — \  2 

x  +y  =  D 
z  =  eH. 


(4.1) 

(4.2) 


where  e  =  =bl.  Equation  4.2  yields 


t  =  {eH  -  az)/{voz  -  viz). 


Equation  4.1  rewrites  to 


(ax  +  t(kvox  -  vix ))2  +  (t(kvoy  -  viy ))2  =  D2 , 


i.e., 


^  i^ox  d"  “1“  [2(Ucc  tTix)tv0x  V0yViy\k  +  [((2#  ^  V^y  D  ]  —  0. 

We  solve  this  equation  for  &. 

2.  Points  on  the  lines  We  have 


x  =  D2 /ax 

V  =  -eDy/al  -  D2/ax, 


(4.3) 

(4.4) 


where  e  —  ±1.  From  Equation  4.3, 

&x  H-  t(kvox  Vix}  —  kD  / (2^ 

fa  —  (D  CLX)  /  {tCLxVox)  +  Vix/vox. 


Equation  4.4  yields  to 

(-D  dx}v0y  j  ( CLxVox )  +  tXfa, / Vqx  's/ 6^  D ^  / (2a, .  (4.5) 
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Fig.  3.7.  Target  set.  Case  az  =  —H  and  —ax  >  D 


Fig.  3.8.  Target  set.  Case  az  =  H  and  —ax  >  D 
We  solve  Equation  4.5  in  t: 

t  =  {-£Dvoxy/a%  -  -D2  +  (a2x  -  D2)voy) / {ax(vixvoy  -voxviy)). 

Thus,  we  deduce 

k  =  {viyC^x  -  -D2  -  vix£D)/(v0yCa?x  -  £>2  -  voxeD). 

We  analyze  the  singularities. 

•  Case  voz  =  Viz.  The  vector  vQ  —  v\  is  horizontal  and  will  remain  horizontal  if  we  change  the 
ground  speed  of  the  ownship.  Solutions,  if  any,  are  tangential  to  the  circle  and  also  belong  to 
P'2  and  those  are  handled  in  case  2. 

Therefore,  there  are  no  solutions. 

•  Case  eDvox  =  voy^Ja2x  —  D 2.  No  solution,  since  t  —  0. 

•  Case  —  ax  <  D.  No  solution,  since  the  set  P'2  is  empty. 

•  Case  —  ax  =  D.  There  is  a  single  line  and  the  relative  initial  position  of  the  ownship  is  on 
that  line.  The  only  way  to  reach  a  target  point  on  that  line  is  if  v'  is  vertical  or  null.  If 
VixV0y  ~  voxViy  =  0  then  there  is  a  solution:  v'ox  =  Vix,v'oy  =  v\y.  Otherwise,  there  is  none. 

•  Case  Vixvoy  —  voxViy  =  0.  The  horizontal  components  of  v\  and  v*Q  are  parallel,  these  two 
horizontal  components  cannot  be  the  same  as  the  aircraft  are  predicted  to  be  in  conflict  and 
— ax  >  D ,  thus  the  relative  speed  v  cannot  be  vertical  or  null.  As  the  horizontal  part  of  vQ  is 
not  0,  the  only  way  to  change  the  direction  of  the  relative  speed  vector  by  changing  the  ownship 
ground  speed  is  to  take  v'ox  =  Vix,v'oy  =  V{y.  This  way,  v'  is  vertical  or  null.  If  voz  ^  Viz  then 
v'  is  vertical  and  since  — ax  >  D,  this  is  not  an  optimal  solution.  If  voz  =  Viz  then  we  have  the 
solution  v'  =  0. 

•  Case  vox  —  0.  Since  Vixvoy  —  voxViy  ^  0,  we  have  V{x  ^  0  and  voy  ^  0.  We  also  have 

X  —  tJLx  tVix 

and  x  reaches  D2 /ax  at  t  —  ( a 2  =  D2) / (axVix)  independently  of  k.  At  that  time,  we  must  have 

y  =  -eDsJa?x  -  D2/ax 
(ax  ~  D2)/(axvix)(kvoy  -  viy)  =  -eD<Jal  -  D2/ax 

k  —  (v iy  £DVix  j  {voy ~D^ ) ) . 
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Fig.  4.1.  Superset  of  target  points 


This  formula  is  a  particular  case  of  the  general  one. 

3.  Special  case  v'  =  0.  The  only  way  to  reach  v'  =  0  by  changing  the  ownship  ground  speed  is  to  have 
VixV0y  ~  voxViy  =  0  and  voz  =  ViZ.  This  is  already  included  as  a  limit  case  above. 

4.2.  Ground  Track  Change  Only.  Since  ground  speed  and  vertical  speed  are  constant,  we  have 

vox  +  voy  =  vlx  +  vly  (4-6) 

v'oz  =  voz  (4.7) 

and  we  must  determine  the  possible  v'ox  and  v'oy  such  that  at  some  time  £,  the  point  (x,  y ,  z)  =  ( ax  +  t(v'ox  — 
vix),t(Ky  -  viy),az  +  t(Kz  -  viz ))  ^  on  the  target  set. 

1.  Points  on  the  circles  P[.  We  have 

x2  +  y2  =  D 2  (4.8) 

*  =  sH.  (4.9) 


Equation  4.9  gives  t  =  ( sH  —  az)/(yoz  —  Viz).  Equation  4.8  rewrites  to 

2t2voyviy  =  («*  -  tvix )2  +  t2v2iy  +  t2v'2x  +  t2v'2y  -  D2  +  2 (ax 
=  (ax  -  tvix)2  +  t2v2y  +  t2v2ox  +  t2v2oy  -  D2  +  2 (ax 
—  E  2  ((I  ,,  tvj  ,. )  tv OT . 

where  E  =  (ax  —  tvix )2  +  t2v2y  +  t2v2ox  +  -  D2.  Therefore, 

^v'oyvly  =  (E  +  2(ax  -  tvix)tv'ox)2 . 

Using  Equation  4.6,  we  get 

^v2v{v2ox  +  v2oy  -  v'2x)  =  (E  +  2 (ax  -  tvix)tv'ox)2 

Hence, 

4t2((a„  -  tvix)2  +  t2v2y)vl2x  +  4 E(ax  -  tvix)tv'ox  +  E2  -  4 t4v2y(v2x  +  v2y)  =  0.  (4.10) 


-  tvix)tv'ox 
~  tvix)tv'ox 
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To  solve  Equation  4.10  on  v'ox,  we  compute  its  discriminant  A  =  lQv^t^A',  where 

A'  =  -E2  +  4 t2(v2ox  +  v%){[ax  -  tvlx)2  +  v2yt2). 

If  A'  >  0,  the  solutions  are 

v'ox  =  (~E(ax  ~  tvix)  +  s'viytVK')/ (2 t((ax  -  tvix )2  +  t2v2y)),  (4.11) 

where  s'  —  d=l.  Then,  we  get 

voy  —  (E  +  2 (ax  —  tViX)tvox) / (2 1  Viy)  (4-12) 

=  ( Etviy  +  e'(ax  -  tvix)\/A')/(2t((ax  -  tvix )2  +  t2v2iy)).  (4.13) 

We  analyze  the  singularities. 

•  Case  voz  =  ViZ.  The  vector  vQ  —  vl  is  horizontal  and  will  remain  horizontal  if  we  change  the 
ground  track  of  the  ownship.  Solutions,  if  any,  are  tangential  to  the  circle  and  also  belong  to 
P2  and  those  are  handled  in  case  2. 

•  Case  =  0  and  ax  ^  tvix.  The  two  solutions  of  the  equation  for  v'ox  given  by  the  formula  4.11. 
There  are  two  solutions  for  v'oy 

^ oy  ^  \J ^ox  "b  ^ oy  ^ ox  ' 

In  fact,  this  formula  is  a  particular  case  of  the  general  one. 

•  Case  =  0  and  ax  =  tvix.  At  time  £,  the  intruder  will  be  where  the  ownship  is  at  time  0. 
Changing  the  ownship  ground  track  does  not  affect  the  horizontal  distance  between  the  aircraft 
at  time  t  (that  will  be  t^v^x  +  v2oy  in  all  cases).  Therefore,  it  does  not  help  to  solve  the  conflict. 

2.  Points  on  the  lines  P'2.  We  have 

x  =  D2/ax  (4.14) 

y  =  -eDs/al-D*lax.  (4.15) 

If  —  ax  7^  D  then  v'ox  =  Vix  is  not  a  solution  and  Equation  4.14  gives  t  in  function  of  v'ox 

t  D2-4 
ax(Kx  -vixy 

Equation  4.15  rewrites  to 

(D2  ~  a2x)l{ax{v'ox  -  vix))(v'oy  -viy)  =  -eD^a2x  -D2/ax 
( V'oy  ~  Viy)/(v'ox  ~  Vix)  =  eD / yj a2  -  D2 

v'oy  =  e£>/\Al  -  D2{v'ox  -  Vix)  +  Viy 

Replacing  v'oy  in  Equation  4.6,  we  get 

(eD/yJa2  -D2(v'ox  -  vix)  +  viy)2  +  vl2x  =  v20X  +  v2oy.  (4.16) 

We  solve  Equation  4.16  and  get  v'ox  and  then  v'oy. 

We  analyze  the  singularities. 

•  Case  —  ax  <  D.  No  solution  since  the  set  P'2  is  empty. 

•  Case  —  ax  =  D.  The  only  solution  is  v'ox  =  ViX ,  v'oy  =  v\y.  In  this  case,  we  must  have 

vix  +  viy  ~  vlx  +  vly  Hence,  v'  is  vertical  or  0. 

3.  Special  case  v'  =  0.  The  only  way  to  get  v'  =  0  by  changing  the  ownship  ground  track  is  to  have 

vix  +  Viy  —  Vox  +  Voy  an4  voz  =  Viz-  In  this  Case,  we  take  v' o  =  Vi. 
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4.3.  Vertical  Speed  Change  Only.  If  the  target  point  is  on  the  line  P it  is  also  in  P[.  We  only 
need  to  consider  the  cases  where  the  point  is  in  P[  and  v'  =  0. 

1.  Points  on  circles  P[.  Thus,  we  have 

x2+y2  =  D2  (4.17) 

z  =  sH.  (4.18) 

Equation  4.17  rewrites  to 

(ax  +  t{vox  -  vix))2  +  ( t(voy  -  viy))2  =  D2 .  (4.19) 

We  solve  Equation  4.19  to  get  t.  Equation  4.18  yields  to 

v'oz  =  viz  +  ( eH  -  az)/t. 

If  vox  —  Vix,voy  =  Equation  4.19  has  a  solution  only  when  —  ax  =  D  and  in  this  case  there  is  no 
conflict. 

2.  Special  case  v'  =  0.  The  only  way  to  reach  v'  =  0  by  changing  the  ownship  vertical  speed  is  to  have 
Vox  =  vix  and  voy  =  viy.  We  take  v'oz  =  viz. 

5.  A  Prototype  Implementation.  We  have  experimented  this  algorithm  with  a  prototype  imple¬ 
mentation.  The  prototype  is  about  a  couple  of  hundred  lines  of  Java,  containing  only  assignments  and 
conditionals.  The  functions  used  are  the  four  operations  and  square  root,  but  no  trigonometric  functions 
(except  in  the  interface,  to  print  the  ground  track  of  the  aircraft  from  the  computed  Cartesian  coordinates 
of  the  speed  vector).  The  implementation  is  available  at  http : //www.  icase .  edu/~munoz/sources  .html. 

Here  is  a  typical  execution:  We  have  two  aircraft  flying  at  the  same  altitude  with  a  horizontal  separation 
of  10  nm.  In  the  coordinate  system  where  the  intruder  is  at  the  origin  and  the  ownship  at  coordinates 
(—10,0,0),  the  ownship  ground  track  is  0  and  the  intruder  ground  track  is  180°.  The  ground  speed  of  the 
ownship  is  400  nm/h  and  that  of  the  intruder  is  300  nm/h.  The  ownship  is  climbing  at  a  vertical  speed  of 
1000  ft/mn  and  the  intruder  is  descending  at  a  vertical  speed  of  —1000  ft/mn.  The  input  to  the  algorithm 
is  a  file  containing  the  following  information. 

Ground  distance  =  10  nm  Vertical  distance  =  0  ft 
Ownship:  0  deg  400  nm/h  1000  ft/mn 

Intruder:  180  deg  300  nm/h  -1000  ft/mn 

The  programs  detects  a  conflict  and  proposes  five  solutions: 

Conflict  in  the  time  interval  (25.7143,29.1456) 

There  are  5  solutions. 

Modify  GROUND  SPEED  317.5889  nm/h  (TOP) 

Modify  GROUND  TRACK  29.1888  deg  (TOP) 

Modify  GROUND  TRACK  -29.1888  deg  (TOP) 

Modify  VERTICAL  SPEED  1266.8799  ft/mn  (TOP) 

Modify  VERTICAL  SPEED  -3266.8799  ft/mn  (BOTTOM) 

The  first  solution  is  to  reduce  ground  speed  to  317  nm/h.  The  second  and  third  modify  ground  track. 
The  last  two  solutions  modify  the  vertical  speed.  On  the  other  hand,  in  the  first  four  solutions,  the  target 
points  are  on  the  top  circle  of  the  target  set.  In  the  last  solution,  the  target  point  is  on  the  bottom. 
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Notice  that  some  solutions  may  not  be  physically  possible.  For  instance,  the  last  solution  proposes  an 
absolute  change  of  vertical  speed  of  more  than  4000  ft/mn.  The  algorithm  does  not  distinguish  between 
the  solutions.  In  fact,  it  is  intended  to  be  used  in  a  more  general  system,  where  the  choice  of  one  solution, 
among  the  multiple  that  have  been  proposed,  may  use  other  kind  of  information  such  as  type  of  aircraft, 
weather  conditions,  other  potential  intruders,  intent  information,  etc. 

6.  Conclusion.  We  have  given  a  complete  and  rigorous  analysis  of  tactical  detection  and  resolution  of 
air  traffic  conflicts  in  the  3-dimensional  space  and  described  a  new  CD&R  algorithm  that  produces  a  set  of 
solutions.  Each  solution  proposed  by  the  algorithm  is  a  constrained  single  maneuver  that,  when  performed 
by  the  ownship,  solves  the  conflict  without  collaboration  of  the  intruder  aircraft.  Experiments  have  indicated 
that  our  algorithm  always  yields  at  least  two  solutions.  After  thousands  of  randomly  generated  examples 
the  average  was  three  solutions  per  conflict. 

Although  the  algorithm  only  uses  state  information,  it  can  be  integrated  within  a  more  general  system, 
such  as  AOP,  to  detect  and  solve  conflicts  in  piecewise  linear  flight  plans.  It  is  well  suited  to  serve  this 
purpose.  First,  it  is  efficient.  Particularly,  it  does  not  contain  loops  nor  calls  to  trigonometric  functions. 
Moreover,  intent  information  can  be  used  to  chose  among  the  multiple  solutions  that  are  proposed.  We  plan 
to  pursue  this  direction  of  research  in  future  work. 

In  the  near  future,  we  will  formally  verify  the  correctness  of  the  algorithm  in  the  PVS  [10]  seek  to  prove 
that  the  solutions  proposed  indeed  solve  the  conflict  and  that  there  is  always  at  least  one  solution  proposed, 
and  extend  our  approach  to  deal  with  multiple  aircraft. 

Acknowledgment.  The  authors  are  grateful  to  Rick  Butler  for  his  helpful  comments  on  preliminary 
versions  of  this  manuscript. 
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